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<N : INTRODUCTION 

^ I The wealth and quality of new experimental findings and the theoretical progress in the 

^ ■ construction and use of Effective Field Theories (EFTs) of QCD are among the reasons 

^ '. of the heavy quarkonium physics renaissance witnessed during the last years. In order 

I to keep track and make immediately available to a larger community the progress in the 

^ ■ field, experimental and theoretical physicists have gathered in the last three years to form 

, a Quarkonium Working Group [1]. The offspring of this collaboration have been three 

O I workshops, a school and a newly issued CERN Yellow Report 1121]. In the following we 

^ ' will review some recent results and open problems in heavy quarkonium physics. We 
refer the reader to L2J for exhaustive presentations. 

Oh! 
(D . 

^. SPECTROSCOPY 

> . 
• 1—1 , 

^ . There has been in recent years a renewed interest and a noteworthy progress in the 

^ i calculation of the heavy quarkonium spectrum from perturbative QCD. The progress 

■ - - ' comes essentially from three directions: 

(1) progress in the construction of EFTs of QCD suitable to describe non-relativistic 
bound state systems. These have helped to organize higher-order calculations and to 
factorize high-energy perturbative and low-energy non-perturbative contributions. 
For a recent review we refer to |0]. 

(2) fixed order calculations f4','5',10,[ill[iS[iil[il[il. 

(3) resummation of large contributions (large logs or large contributions associated to 
renormalon singularities) [14. 15. 16. 17.18.191 . 

By definition, a perturbative treatment of the bound state is possible if the momentum 
transfer scale p of the heavy quarks in the bound state is much larger than Aqcd- Two 
situations may occur under this condition boll . Let us call E the typical kinetic energy of 
the heavy quark and antiquark in the centre-of-mass reference frame: in a non-relativistic 
bound state p ~ mv ^ £ ~ mv^, v being the heavy-quark velocity in the bound state. 



The first situation corresponds to quarkonium states for which E > Aqcd- Under this 
circumstance the heavy-quarkonium potential is purely perturbative; we may call this 
case Coulombic. In the case E ^ Aqcd non-perturbative contributions are encoded 
into local condensates, in the case E ~ Aqcd into non-local ones (see also llllll and 
references therein). The second situation corresponds to quarkonium states for which 
p ^ Aqcd ^ E. Under this circumstance the potential contains a perturbative part and 
short-range non-perturbative contributions. We may call quasi Coulombic the case in 
which non-perturbative contributions to the potential turn out to be small and may be 
treated perturbatively in the calculation of observables. 

Clearly, perturbative calculations applied to heavy-quarkonium ground states are on 
a more solid ground that those to higher resonances, and those applied to bottomonium 
more than those to charmonium, since the (quasi-)Coulombic case is more likely to be 
realized. Perturbative determinations of the T(15) and 7/i// masses have been used to 
extract the b and c masses (see Masses). The main uncertainty in these determinations 
comes from non-perturbative contributions. Information on their size can come from 
other ground-state observables, like the hyperfine splittings (see Hyperfine splittings), 
and the mass (see Be mass). Recent studies of higher quarkonium states in perturba- 
tion theory may be found in |4, ,22j|,|23||,|5||,|24j,i25»il2|]. Likely only some of the lowest 
bottomonium resonances may be treated consistently as Coulombic or quasi-Coulombic 
bound states. It is surprising, therefore, that some of the gross features of the bottomo- 
nium spectrum, like the equal spacing of the radial excitations, may be qualitatively 
reproduced by pure perturbative calculations (see Higher bottomonium states). 



TABLE 1. Collection of recently obtained values of 
mp(mp) and mf^(mp) from the T(1S) and 7/^ masses. 



Ref. 


order 


mf^imf^) (GeV) 


[27J 

[25] 
[26] 
[30] 
[31] 


NNLO 
NNLO 
NNLO 
NNLO 
NNNLO 
NNNLO 
NNNLO 


4.24 ±0.09 
4.21 ±0.09 
4.210 ±0.090 ±0.025 
4. 190 ±0.020 ±0.025 
4.349 ±0.070 

4.20 ±0.04 
4.241 ±0.070 




Ref. 


order 




\M 


NNLO 


1.24 ±0.020 



Masses 

Table[T]shows some recent determinations of the b and c masses from the quarkonium 
ground-state masses. Finite charm-mass effects have been included in the analyses of 
[25, 30, 31]. With the exception of |27] the conversion from the threshold (or pole) 
masses to the IMS masses has been performed at three-loop accuracy. All NNNLO 
analyses are only partially complete, since the full three-loop static potential has not been 
calculated yet. Different schemes have been used to implement the leading renormalon 



cancellation. This may explain some of the differences between the results. The most 
substantial discrepancy, which is between the result of |26J and all the other ones, 
may be possibly ascribed to the use of the on-shell scheme, as well as to the way 
non-perturbative effects have been implemented. In i24l1 the charm MS mass has been 
extracted from the J /m mass. 

In the above analyses non-perturbative effects constitute the major source of uncer- 
tainty together with possible effects due to subleading renormalons. Indeed, in the situa- 
tion E ~ Aqcd non-perturbative corrections to the spectrum (~ Aq^^/ ~ mv^) and the 

subleading renormalon in the mass (~ J^Q^/m ~ mv^) are parametrically of the same 
order as NNLO corrections in the perturbative expansion. As long as non-perturbative 
effects will not be incorporated in a quantitative manner, perturbative calculations be- 
yond NNLO will not improve the determination of the heavy-quark masses from the 
quarkonium system. 

The dominant non-perturbative correction to the mass of a quarkonium state \n) of 
LO binding energy E^"^ is encoded in the expression 

dE„ = -i^ J?(n|re"(^« -''"^rln) (E(0 E(0))(Ai), 
9 Jo 

ho = p-^/m + as/ (6r) being the LO octet Hamiltonian and /i a factorization scale lIsLlloll. 
A precise determination of this formula would need having an accurate determination of 
the chromoelectric correlator. In this respect the available data are puzzling. Quenched 
lattice calculations of gluelump masses indicate an inverse correlation length for the 
correlator (E(?) E(0)) of about 1.25 GeV [32, 33]. This scale is not only larger than the 
scale E but also of the same order as the momentum-transfer scale p, if p is identified 
with mfe/2 X 4/3 X as{p) (e.g. mt, ~ 4.7 GeV and a^ip) ~ 0.4). This may potentially 
question the whole perturbative approach to the ground state of bottomonium. On the 
other hand this approach is supported by the fact that it provides a value of the b mass 
consistent with those obtained from low moments sum rule calculations or from the B 
system ll34ll . Only a better determination of the chromoelectric correlator can solve the 
conundrum. A possible solution could come, for instance, from 

• unquenching, if it lowers the value of the inverse correlation length as calculations 
with cooling techniques seem to indicate [35]. 

• a particular parametrization of the chromoelectric correlator that sets in the non- 
perturbative behaviour at a scale ji lower than the inverse correlation length. 

The assumption on the Coulombic behaviour of ground- state quarkonia may be also 
tested on other observables, like the Be mass and the hyperfine splittings. 



Be mass 

Table |21 shows some recent determinations of the Be mass in perturbation theory at 
NNLO accuracy compared with the value of the Be mass in f36] and in a very recent 
lattice study 1.37.1 . The b and c pole masses are expanded in terms of the T(15) and 



7/i/A masses respectively in fss"], and in terms of the MS masses in ll24l Esll . In 1I25I1 
finite charm-mass effects are included. Since the charm quark effectively decouples in 
the scheme of Ref. [38], we may consider the results of Refs. IbsIi and [25 ] at the same 
level of theoretical accuracy and the difference between them as due to higher-order 
corrections beyond NNLO accuracy. 



TABLE 2. Different perturbative determinations of the mass compared with the experi- 
mental value and a recent lattice determination. 



Be mass (MeV) 


State 


m (expt) 1 


[37] (lattice) [38] (NNLO) 


[24] (NNLO) 


[25] (NNLO) 




6400(400) 


6304 ±121^2 6326(29) 


6324(22) 


6307(17) 



In Fermilab Today, December 3, 2004, the CDF collaboration has announced the 
preliminary results of a search for the Be, using decays into a.J/\j/ and a charged pion. 
Their preliminary result is Mg^, = 6287 ± 5 MeV. This value, if confirmed, would support 
the assumption that non-perturbative contributions to the quarkonium ground state are of 
the same magnitude as NNLO or even NNNLO corrections, which would be consistent 
with aE > Aqcd power counting. 

Hyperfine splittings 

Charmonium and bottomonium ground state hyperfine splittings have been recently 
calculated at NLL in [Ts']. The result in the charmonium case is consistent with the 
experimental value. This, again, supports the assumption that non-perturbative contribu- 
tions for the quarkonium ground state are consistent with a E > Aqcd power counting. 
In the bottomonium case there are not yet experimental data. 

We note, however, that NNLO corrections may be potentially important for this 
observable, also considering that at this order new type of corrections will appear for 
the first time: 

• non-perturbative corrections, which are parametrically of NNLO in the situation 
Aqcd ~ E. 

• corrections to the mass. In order to show the possible impact of this type of 
corrections, in Fig.[T]we display the hyperfine splittings at NLO accuracy calculated 
using 1/2 of the bottomonium and charmonium vector ground state masses and 
the MS masses. The difference between the two determinations (taken around 
the maximum in ju) is about 70% in the charmonium case and about 30% in the 
bottomonium case. 



Higher bottomonium states 

Higher bottomonium resonances have been investigated in the framework of pertur- 
bative QCD most recently in [24, "25113^. The calculation of [24] is accurate at NNLO, 
that one of [25] includes finite charm-mass effects, while [i39i] is a numerical calculation 




1.5 2 2.5 3 3.5 



IX (GeV) IX (GeV) 

FIGURE 1. Hyperfine splittings of the IS charmonium (left) and bottomonium (right) states at NLO 
in the 15 mass scheme, nic = nij, — My(j5)/2, (lower curves) and in the MS mass scheme, 

ruc = mp(m^)(l +4a,/{3n)), nit = m^(mp)(l + 4as/(37r)), (upper curves) as a function of the 
normalization scale pi. Os runs at 4-loop accuracy with 3 massless flavours. Left plot; m^^{m^^) — 1.28 
GeV, as(My/^/2) 0.352, as(mp(m^)) = 0.405. Right plot: mf^{mf^) = 4.22 GeV, as(Mx(is)/2) = 
0.209, as(mfS(,„MS)) = 0.218. 

TABLE 3. Masses of bb states. The result of Ref. i24ll comes from a full perturbative 
calculation up to 0{ma^) without finite charm-mass corrections; the result of Ref. f25| 
comes from a full perturbative calculation up to 0{ma^) including finite charm-mass 
corrections; the result of isoil incorporates full corrections up to 0{m af) in the individual 
levels and full corrections up to 0{ma^) in the fine splittings, includes finite charm-mass 
corrections and also depends on some assumptions about the long-range non-perturbative 
behaviour of the static potential. Numbers without errors are those without explicit or 
reliable error estimates in the corresponding works. 

bb states 



State 


expt 


[24J 


[251 


[39] 




9460 


9460 


9460 


9460 




9913 


9916(59) 


10012(89) 


9956 




9893 


9904(67) 


10004(86) 


9938 




9860 


9905(56) 


9995(83) 


9915 




10023 


9966(68) 


10084(102) 


10032 




10269 




10578(258) 


10270 




10255 




10564(247) 


10260 


23Po 


10232 


10268 


10548(239) 


10246 


3^51 


10355 


10327(208) 


10645(298) 


10315 



of the spectrum that includes also NLO spin-dependent potentials. This last calculation 
also relies on some assumptions about the long-range behaviour of the static potential. 
The results are summarized in Table 

The surprising result of these studies is that some gross features of the lowest part of 





H (GeV) 



FIGURE 2. E 
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■"Xbii^P) (^^S^^) versus the normalization scale jj.. 



The light band shows the LO expectation, the dark one the NLO one. The widths of the bands account 
for the uncertainty in as (Mz) = 0. 1 1 87 ± 0.002 1 36] . The horizontal bands show the experimental values 
is^.FromyjI. 



the bottomonium spectrum, like the approximate equal spacing of the radial levels, is 
reproduced by a perturbative calculation that implements the leading-order renormalon 
cancellation. If this is coincidental or reflects the (quasi-)Coulombic nature of the states 
will be decided by further studies. A recent NLO calculation of the IP bottomonium 
fine splittings has been performed in [Ts]. The results are summarized in Fig. El Figure 
01 plots p, the ratio of the fine splittings considered in Fig. El as a function of the 
normalization scale jU (see i40] and references therein for early studies of this quantity). 
It seems to indicate either the existence of large NLL/NNLO corrections (as it happens 
in the hyperfine splittings of the 15 levels) or sizeable non-perturbative corrections, 
somehow hidden in the error bands of Fig.O 
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FIGURE 3. p = (£^62 ^ ^Xbi) li^Xbi ^^Xbo) versus the normalization scale jj.. The horizontal line at 
0.8 corresponds to the LO expectation, the curve to the NLO one. The band accounts for the uncertainty 
in (Xs{Mz) = 0. 1 187 ± 0.002 1 36]. The horizontal band shows the experimental value From O. 
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FIGURE 4. Different determinations of as{Mz) (left) and of as(/i) at different scales jj. (right). In both 
figures the determinations outside the band come from Y decay data. From i3^ . 



DECAY 

T inclusive decays 

The NRQCD factorization fomiulas state that inclusive decay widths to light hadrons 
(l.h.) may be written as sums of products of matrix elements and imaginary parts of 
Wilson coefficients of 4-fermion operators ll4lll . In the case of the T system they read: 



r{T^ l.h. 



Im/i(35i)(Oi(35i))T 



+lmfses,){OsCSo)U, + ^^^^^^^{OseS,)) 



Vh 



+ ■ 



+ ■ 



The use of these formulas (and similar ones for the electromagnetic decay widths) to 
extract (Zg at the Z^-mass scale has turned out problematic, as shown by Fig.|4l It seems 
that uncertainties have been underestimated [2]. The reasons may be the following: 

• poor knowledge of the matrix elements. Note that usually octet matrix elements 
have been neglected in the analyses. This may not be justified considering that 
Im/i(^5i) and Imgi(-^5i) are suppressed by an extra ttg with respect to most 
of the Im/g coefficients (see [421 and references therein) that may compensate, 
under some circumstances, the suppression in the velocity expansion of the octet 
operator matrix elements. The poor knowledge of the matrix element (^i) has 
been indicated as one of the major source of uncertainty in lG^3l1 . Studies of decay 



matrix elements on the lattice can be found in 147. '45'. ^4^. Factorization formulas 
for decay matrix elements in pNRQCD have been derived in ll47..i48i1 . 
• large higher-order corrections to the Wilson coefficients ll49l] . 

TABLE 4. Comparison of ratios of %c] partial widths. The experimental values PDG 2004 and 
PDG 2000 are obtained from the world averages of |36] and |50] respectively. The chosen ratios 
do not depend at leading order in the velocity expansion on octet or singlet matrix elements. 
The LO and NLO columns refer to a leading and next-to-leading order calculation done at the 
renormalization scale 2mc with the following choice of parameters: nic = 1.5 GeV and a^iTiric) = 
0.245. From |2]. 



Ratio II [36] (PDG 2004) 


[50] (PDG 2000) II LO 


MT O 

iNJ_,W 


r(z.o ^ 77) 
r(z<2 ^ 77) 


5.1±1.1 


13±10 


« 3.75 


«5.43 




410±100 


270±200 


w347 


« 383 


r(z.o ^ 77) 


r(Zco^/./i.)-r(z.i^/A) 


3600±700 


3500±2500 


w 1300 


« 2781 


r(ZcO 77) 


r(z,o^/./i.)-r(x,2^/./2.) 


7.9±1.5 


12.1±3.2 


«2.75 


«6.63 




nXc0^l.h.)-T{Xci^l-h.) 


8.9±1.1 


13.1±3.3 


«3.75 


«7.63 


nXc2^l-h.)-T{Xci^l.h.) 



Xc inclusive decays 



In the case of Xc inclusive decay widths the NRQCD factorization formulas read: 

2 



TiXcj-^l.h.) = 9lmfiCPj) 



ntrr 



+ 



+ ■ 



7 = 0,2, 



where ^^^^(0) is the derivative of the radial part of the Xcj wave function at the origin. 
The displayed terms are accurate at leading order in the velocity expansion. The Wilson 
coefficients are known at NLO accuracy in a^. 



The experimental determination of the Xc decay widths has dramatically improved 
in the last four years mainly due to the measurements of the E835 experiment at the 
Fermilab Antiproton Accumulator. A way to see the impact of these measurements is 
provided by Tab. S where we compare the PDG 2000 Q with the PDG 2004 S 
determinations of different ratios of hadronic and electromagnetic widths. There have 
been sizable shifts in some central values and considerable reductions in the errors. In 
particular, the error on the ratio of the electromagnetic ;fcO and Xci widths has been 
reduced by about a factor 10, while in all other ratios the errors have been reduced by 
a factor 2 or 3. The considered ratios of hadronic and electromagnetic widths do not 
depend at leading order in the velocity expansion on any non-perturbative parameter. 
Therefore, they can be determined in perturbation theory only. The last two columns 
of Tab. m show the result of a LO and NLO calculation respectively. Clearly the data 
have become sensitive to NLO corrections and may be used, in principle, to determine 
tts at the charm-mass scale. Before this, a necessary step is the calculation of the decay 
widths at next-to-leading order in the velocity expansion, since these contributions are 
potentially of the same magnitude as NLO corrections in the Wilson coefficients (see 
Ref. llsill for a calculation in the electromagnetic case). 

PRODUCTION 

In this section we briefly mention two of the main open problems in our understanding 
of charmonium production. We refer to L2J for a proper treatment. 



Charmonium polarization 

Quarkonium production at large pj is dominated by gluon fragmentation. NRQCD 
predicts that the dominant gluon-fragmentation process is gluon fragmentation into a 
quark-antiquark pair in a color octet state. At large pj the fragmenting gluon is 
transversely polarized. In the standard NRQCD power counting it is expected that the 
octet quark-antiquark pair keeps the transverse polarization as it evolves into a 5-wave 
quarkonium. Different power countings [52] or higher-order corrections may somehow 
dilute the polarization, which, however, is expected to show up in high pj data. The 
present Tevatron data [53] do not seem to confirm this expectation. The uncertainties 
are, however, too large to make any definite claim. The issue here is mainly experimental 
and more precise Tevatron data are eagerly awaited. 



Double charmonium production 

The most challenging open problem in charmonium production concerns double 
charmonium production in e+e^ annihilations. The Belle data of inclusive 7/1// + cc 



production, 

C7(e+e" — >7/V/ + cc) 



0.82±0.15±0.14, Belle [54], 



a{e+e- -^J/\j/+X) 
and of exclusive J/xir + rjc production, 

a(e+e" ->7/V/+77c)Br(ccres ^> 2 ch. track) =25.6 ±2.8 ±3.4 fb, Belle JH], 
are far by almost an order of magnitude from the available theoretical predictions: 

a{e+e- ^J/\j/ + r]c) = 2.31 ± 1.09 fb, 

On the other hand, the present upper bound of Belle on 7/i//±7/i// production is 
consistent with theoretical expectations: 

a(e+e" ^7/v/ + 7/v/)Br(ccres ^> 2 ch. track) < 9.1 fb. Belle S, 



(j{e+e- ^J/\j/ + J/\j/) = 8.70±2.94fb, 116 iL 16211. 



Even if independent data by BaBar would be most welcome, the issue here seems to be 
mainly theoretical and the data may be the signal of some new production mechanism. 
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